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I. ABSTRACT 

Theoretical models of an axisymmetric laminar diffusion flame are 
discussed, with an emphasis on the behavior of such flames at increasing 
pressures. The flame-sheet or Burke-Schumann model (in terms of Bessel 
functions) and various boundary- layer numerical solutions are presented 
and their results compared with experimental data. The most promising 
theoretical model combines the numerical flow field solution of the 
Patankar-Spalding computer code with the Pratt-Wormeck chemical reaction 
subroutine. The flame shapes for pressures of 1 , 5, 10, 20, and 50 
atmospheres have been computed and agree remarkably well with exper- 
imental data. There is a noticeable shape change with pressure, believed 
to be a result of buoyancy effects. The chemical concentration profiles 
do not exhibit much dependence on pressure, a reflection of the fac^ 
that only one chemical mechanism was utilized at all pressures. 
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II. INTRODUCTION 

The problem of pollutant formation in turbine combustors is currently 
receiving a great deal of warranted attention, especially the pollution 
associated with aircraft turbine engines in and around airports. Of the 
various pollutants discharged from these engines, two of the most harmful 
are nitric oxide and nitrogen dioxide, commonly referred to as NO . 

Since actual combustors may operate at pressures up to 30 atm, it is of 
special interest to consider the effect of pressure on the formation of 
NOj^ during the combustion process, which has been done experimentally by 
Miller and Maahs [1] at NASA-Langley. 

In their experimental investigation, an intentionally simplified 
combustion system— a laminar diffusion flame— was chosen to model the 
actual combustion process. The diffusion, or non-premixed, flame was 
selected over the premixed flame because it is an experimentally simpler 
system and because larger NO concentrations are expected with the 
diffusion flame. The laminar flow regime was selected in order to 
eliminate fluid flow effects as much as possible so as to concentrate on 
chemical effects. 

The structure of laminar diffusion flames has been the subject of 
investigation for a great many years. The first analytical study of 
diffusion flames was made by Burke and Schumann [2] in 1928. They 
introduced the assumption that the zone of combustion might be approxi- 
mated as a surface, i.e. fnfinitely thin. It would then follow that the 
associated reaction rates would be infinite. Burke and Schumann were 
concerned with a cylindrical laminar fuel jet discharging into an air 
stream moving with the same velocity as the fuel; they proceeded to 
calculate the flame shape for various initial conditions. The same 
assumption has been utilized to predict the length of both laminar and 



turbulent open flames. CHottel and Hawthorne [3], Wohl , Gasley, and 
Knapp [4], Barr [5]). 
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Zeldovich [6] was among the first to consider the structure of a 
reaction zone of finite thickness and hence the effects of finite reaction 
rates. Similar studies have been done by Spalding [7] with an emphasis 
on the phenomenon of extinction. Linan [8] has also studied the effects 
of finite chemical reaction rates by means of a boundary- layer-type 
solution, while many have applied perturbation techniques to the problem 
(Fendell [9], Kassoy and Williams [10], Chung, Fendell, and Holt [11]). 

The computer simulation and prediction of chemically reacting 
flows, an extremely difficult problem because one must describe the 
interaction of many reacting chemical species with fluctuating temperature 
and velocity fields, is nevertheless a method of solution that is being 
used to a greater and greater degree as computer software/hardware 
facilities become more advanced. 

There exist at present two well-defined bodies of literature on the 
numerical (computer) solution to chemically reacting flows: 1) fluid 

flow with infinite-rate chemistry (chemical equilibrium), and 2) finite- 
rate chemically reacting flow without mixing or with infinite-rate 
mixing. Examples of such works are those by Patankar and Spalding [12] 
and Pratt [13], respectively. Pratt and Wormeck [14] have then attempted 
to combine aspects of both bodies of literature to achieve a scheme that 
could model the phenomena of a diffusion flame-a laminar finite-rate 
chemically reacting flow field. 

It is the purpose of this study to construct an analytical model of 
the laminar diffusion flame described by Miller and Maahs [1] and to 
compare experimental and theoretical results concerning the effects of 
high pressure on flame shape, temperature distribution and NO formulation. 
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Sections III and IV describe the physical system and experimental results. 
Section V introduces the governing equations and Section VI describes 
various theoretical models and associated results. 
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III. DESCRIPTION OF THE PHYSICAL SYSTEM 

For purposes of completeness a summarized description of the 
physical system employed at NASA-Langley is included. For a more 
detailed description. Reference [1] should be consulted. 

A diagram of the burner and collection system is shown in Figure 1. 
Methane gas issues from a cylindrical tube into a concentric cylindrical 
tube through which air is flowing. The fuel tube has an inside diameter 

of 3.06 mm and the air tube has an inside diameter of 20.5 mm. Fuel and 
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air flow rates are held constant at 41.8 standard cm /min (seem) and 
2450 seem, respectively. 

In order to obtain a flat velocity profile in the fuel and air 
streams as they exit from their respective ducts, porous disks of 
sintered stainless steel were installed near the duct exits. Ignition 
is accomplished by means of an electrical hot-wire igniter. During 
operation the igniter is positioned 10 to 20 mm above the burner^ once 
ignition is attained the igniter is removed from the flow. 

The combustion products exit from the top of the quartz chimney and 
are swept into a total -sample collector by the nitrogen pressurizing 
gas. The nitrogen flow rate is kept low enough so as not to excessively 
dilute the gases to be analyzed. The products are analyzed by a continuous 
infrared COr, monitor and a chemiluminescent NO-NO analyzer. 

^ X 

The measurement of flame temperature is the most difficult and most 
uncertain of any measurements made. The method of two-color pyrometry 
was employed (with the carbon in the flame as the radient emitter) to 
determine a spatial average temperature near the hottest point in the 
flame. 



Hot Hire Igniter 



Figure 1 


Schematic diagram of burner and collection system 





7 


IV. EXPERIMENTAL RESULTS 

* 

A. Flame Shape 

A definite change in the shape of the diffusion flame occurred as 
pressure was increased from 1 to 50 atm. This is shown in Figure 2, 
taken from Reference [1]. At 1 atm the sides of the flame bow outward, 
but this decreased as the pressure was increased above 1 atm. At about 
10 atm the flame exhibited a slight concavity and the diameter decreased 
considerably. A further decrease in flame diameter occurred as the 
pressure increased, the maximum change associated with a pressure of 
approximately 20 atm. 

An increase in the amount of carbon accompanied the change of flame 
shape and at pressures higher than 20 atm only minor changes occurred. 

The height of the flame did not change drastically with pressure, 
as illustrated in Figure 3. The flames were typically 10 mm in height. 

B. Nitrogen Oxide 

The molar emission index of N0„, Im, is defined as the number of 

A 

moles of NO and NOg produced per mole of fuel burned. Figure 4 shows 
that this quantity increases rapidly as pressure Increases, reaching a 
maximum at about 9 atm. Above this pressure the molar emission index 
decreases slowly and continuously until it returns to roughly its original 
value. 

C. Flame Temperature 

A plot of the optically-measured peak temperature as a function of 
pressure is shown in Figure 5, with data limits shown by vertical lines. 
There is a steady decrease of peak temperature with pressure. 




Ambient pressure , atra 


Figure 2 - Photographs of methane-air diffusion flames for various eunbient pressures. 

Methane flow rate, tl.8 seem; air flow rate, 2t50 seem; burner diameter, 3.06 mm. See Ref. 1 
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V. GOVERNING DIFFERENTIAL EQUATIONS 


The equations of motion reflecting the conservation of mass, momen- 
tum, energy, and mass of each component for a chemically reacting gas 
mixture are: 

Conservation of Mass 

1^ + div p V = 0 

where p = density V = velocity vector 
Newton's Second Law 


= -Vp + div T 

where -^ = substantial derivative 
p = pressure 

T = viscous stress tensor 

Conservation of Energy 

Dh Dp , . . — 

= pR + T : e . div q 

where 


at 


+ V 


h = enthalpy 

e = rate of strain tensor 

q = heat flux vector 

S pecies Conservation 
Onii 

Pd^= - div j. + R. 


V 


origfnal: page is 

OF POOR QUALITY 


where 

m^. = mass fraction of species i == — 

= mass flux vector 

R. = rate of creation of species i by chemical reactions 
In order to close the system the following constitutive relationships 
and equations of state are introduced: 



13 


T “ 2\xe + Xdiv Y n 

where vi + A are the 1st and 2nd coefficient of viscosity and 
n 15 the unit tensor 
e = 1/2(VT + VT^) 

where T* is the mass flux vector and D^^'s are the diffusion coefficients 

1 V 1 J 

q = - 1^77 + 2h-J^ ^ is the heat flux vector T = temperature 


mi 

p = pRT 2 ^ 
i i 


k = coefficients of conductivity, 

h = h, = + J C 


dT 


where = molecular weight of species 1. 


The above equations are non-linear second order partial differential 
equations. The two main types of solutions are analytical and numerical. 

Analytical methods are those that express the solution as algebraic 
formulae, power series, etc. Methodologies such as separation of variables 
series expansions, and coordinate transformations are employed to obtain 
such solutions. These types of solutions are usually limited to linear 
or linearized differential equations. 

Numerical methods are those that express the solution in the form 
of numbers related to the particular circumstances of the flow. The 
differential equations are transformed to algebraic "finite-difference" 
equations and the values of all dependent variables are calculated at 
the nodes of a grid covering the domain of interest. 

The flow to be considered is in general governed by nonlinear 
equations with the rate of chemical production term (R^. ) being highly 
complex. Simplifying assumptions may be made and various solutions 
obtained. A discussion of these solutions and associated assumptions 
follows in Section VI. 
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VI. THEORETICAL ANALYSES 
A* Burke - Schumann Analytical Solution 

If the following assumptions are made: 

1. steady cylindrical flov/ 

2. inviscid flow 

3. negligible pressure and thermal diffusion 

4. Picks Diffusion Law 

5. single step chemical reaction of fuel and oxidant ^ product 

6. parallel and equal gas velocities 

7. constant mass flux 

8. axial diffusion «< radial diffusion 

9. entire chemical reaction occurs at the flame surface 

then the governing equations may be combined into a single linear partial 
differential equation with appropriate boundary conditions. The solution 
is in terms of Bessel functions and a flame shape may be predicted from 
the analysis. Reference [15] presents a more detailed discussion of 
this solution. 

The predicted flame shape agrees fairly well with the experimental 
evidence, although the predicted flame height at 1 atm (3.6 cm) is 
greater than the experimental. Furthermore, the Burke-Schumann model 
predicts absolutely no change of shape with pressure, hence it is evident 
that a more refined theoretical analysis is necessary. 

B. Edelman - Fortune - Heilerstein Numerical Solution 

The effects of inertia, viscosity, diffusion, gravity, and combus- 
tion on an axisymmetric laminar diffusion flame have been investigated 
via a boundary layer type formulation by Edelman, Fortune, and Weilerstein 
[16], First a von Mises transformation is made in going from axisymmetric 
physical coordinates (x, r), where x is the axial distance and r the 
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radial distances to {x, coordinates. The stream function is defined 


by 


J = 




r ^ 

J purdr 


The governing equations are then non-dimensional ized by means of an 


emperical normalization parameter. 

"It is customary practice to do this (non-dimens ionalizati on) 
by normalizing the variables directly in terms of the boundary 
conditions. However . , . the resulting dimensionless groups do 
not truly characterize the process and that is due to the large 
variation in flow properties throughout the domain of interest. 
Accordingly, the current approach involves the introduction of a 
set of characteristic quantities which are based upon some state 
within the flow that more accurately reflects the flame structure." 
[16] 


The equations are then solved numerically by means of an explicit finite 
difference formulation. 

The Burke-Schumann combustion model, or "flame sheet assumption”, 
is that all the heat is released on the suface where the fuel -air ratio 
is stoichiometric. A more accurate combustion description is offered in 
this model in terms of a basic equilibrium theory for the combustion in 
either the fuel-lean or fuel-rich regime. 

For a pressure of one atmosphere this program predicted a flame 
height of .8 cm and a flame shape that agrees remarkably well with 
experimental data. However, the results for pressures greater than one 
atmosphere could not be substantiated either qualitatively or quantitatively 
by experimental data. This was due to the fact that the analysis in 
Ref. [16] utilizes. a normalization parameter obtained experimentally at 
one atmosphere of pressure. This does not affect the original work 
since the problem considered was buoyancy effects of diffusion' flames , 
all at one atmosphere. The analysis, limited as it is to a pressure of 
one atmosphere, does not therefore lend itself to the problem at hand. 



C. Patankar-Spalding Numerical Solution 


Patankar and Spalding [12] likewise utilize a boundary layer formu- 
lation and a von Mises coordinate transformation so as to put the 
conservation equations in the general form 

|i + (a + bw) d 

3x ^ ^ aaj^3o) ‘ 

where 

X and w are the two independent coordinates 

a and b are arbitrary functions of x 

c and d are arbitrary functions of any dependent and independent 

variables 

(f) stands for any one of a set of variables, each possessing its 
own differential equation, e.g., fluid velocity, stagnation 
enthalpy, mass fraction, etc. 

These differential equations are transformed into finite difference 
equations and solved by means of the tri -diagonal matrix algorithm 
(TDMA). 

The chemistry incorporated into this numerical scheme is a very 
simple fuel -oxidant-product equilibrium model where the rate of genera- 
tion of fuel by chemical reaction, R.- , is assumed to obey an "Arrhenius 
type" relation. 

For the case of negligible chemical reactions this program produces 
acceptable results in that the diffusional nature of the axial flow of 
fuel and air is well represented by means of velocity and concentration 
profiles. With the inclusion of limited equilibrium chemistry a flame 
shape can be obtained that agrees qualitatively with that obtained 
experimentally. The theoretical flame height 

at one atmosphere is 8 cm as opposed to an experimental height of .8 cm 
however. Since the chemistry of the flow field under consideration is 
of prime importance, the very simple chemistry incorporated into this 
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model does indeed limit its ability to accurately describe the flame 
under investigation. 

D. Patankar-Spalding and Pratt-Wormeck Numerical Solution 

As stated previously, the hydrodynamic aspects of Spaldings code 
are well -developed and offer a valid model of the diffusional nature of 
the flow. A more realistic chemical model is provided by Pratt and 
Wormeck [14] in their subroutine for the calculation of chemically 
complex equilibrium or non-equilibrium stationary states. 

The scheme employed by Pratt and Wormeck in their determination of 
chemical equilibrium states at prescribed temperature and pressure is 
essentially the same utilized by Gordon and McBride [17]: Gibbs function 

minimization. 

A similiar procedure is employed for the case of non-equilibrium 
(kinetic) stationary states. A chemical-kinetic source term is included 
which takes into account the effects of forward and reverse reaction 
rates, and reaction contact indices, 

A combination of the hydrodynamic code of Spalding and the chemical 
package of Pratt and Wormeck results in a computer code that is capable 
of producing a numerical solution to the hydrodynamic equations of 
motion while taking into account the effects of complex chemical reaction. 
This combination computer code was applied to the problem at hand for the 
two cases of equilibrium and kinetic chemistry. It was found that the 
kinetic solution proved much too costly in computer time (IBM 370) to even 
consider, so the following results are limited to the case of equilibrium 
chemistry. 

The thermochemical data were taken from Gordon and McBride [17]. A 
reaction mechanism consisting of 19 species and 34 reactions was utilized, 
as given in the actual computer listing found in the Appendix. A 



18 


cross "St ream grid of 21 grid points was used. The results were found 
to be fairly sensitive to the input value for the temperature at the 
base of the burner. Since the computer solution is for a steady flow, 
the actual ignition phase is not considered, but a temperature is chosen 
that hopefully represents the physical situation after a steady state 
has been achieved. All the results discussed are for a base temperature 
of 1000°K. 

Figure 6 demonstrates the change in theoretical flame height with 
pressure. As can be seen from a comparison of Figures 6 and 3, this agrees 
very well with the data obtained experimentally. The average computed 
flame height was 8.6 mm, as compared to the value of 10 mm for the 
emperical case. 

Figures 7 and 8 demonstrate the theoretical flame shapes for pressures 
of one and five atmospheres of pressure and ten, twenty, and fifty 
atmospheres, respectively. The flame position at a particular height 
above the burner was assumed to be at the point of maximum temperature for 
that height. 

As can be seen from a comparison of these Figures (7 and 8) 
with the actual flame shapes obtained (Figure 2), there is a remarkable 
agreement between the mathematical model and the emperical data. At one 
atmosphere the computed flame is fairly bulbous in nature, the sides of 
the flame bowing outward away from the burner. As pressure is increased 
the outward bowing of the sides of the flame decreases, and at 10 atmo- 
spheres there is almost no bowing at all. At 20 and 50 atmospheres the 
flame exhibits a concave appearance, which agrees extremely well with 
experimental results. 

This shape change is believed to be a direct result of buoyancy 
effects. As the pressure increases, the difference between the densities 
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of methane and air become much more pronounced, hence the behavior of 
the two now very different gases under the action of gravity, i.e., 
buoyancy, affects the resulting flow field. In addition, the problem at 
hand is one of fairly low velocities. As the pressure is increased 
these velocities become even smaller and the fluid flow more sensitive 
to buoyancy effects. Prior to the inclusion of buoyancy effects in the 
computer program, all flames exhibited a somewhat bulbous shape similar 
to that found at one atmosphere. 

Figure 9 is a plot of approximate maximum flame temperature as a 
function of pressure. Since these values are fairly sensitive to the 
input temperature at the base of the burner, it is difficult to make any 
general statement concerning this data except to say that in comparison 
with Figure 5, the computed flame temperature behaves in qualitatively 
the same manner as the adiabatic flame temperature shown there. 

Figures 10-15 demonstrate various species mass fractions as a 
function of the radial distance outward from the center of the burner. 
Figure 10 is for a pressure of 1 atmosphere at a position of .15 cm 
above the inlet. Figures 11-15 are for a position of 0.5 cm above the 
burner, toward the top of the flame, and for pressures of 1, 5, 10, 20, 
and 50 atmospheres respectively. 

At a vertical distance of .15 cm above the burner there still 

exists some unburned fuel close to the axis of the burner. The maximum 

NO concentration occurs closer to the flame position, as does the 
X 

maximum concentration of COg. 

At a vertical distance of .5 cm above the burner, there is no CH^ 
present and once again the maximum NO production occurs near the flame 

A 
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front. There is essentially no change in the relative magnitude of this 
maximum concentration with pressure, only a change in where it occurs. 
This of course corresponds to the change of flame position with pressure. 


APPROXIMATE MAXIMUM 
FLAME TEMPERATURE^ oK 






SPECIES MASS FRACTION 
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Figure 11 - species mass fraction as function of pressure 
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Figure 13 - species mass fraction as function of pressure 
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VII. CONCLUSIONS 

In order to reflect the physical system under consideration to an 
acceptable degree, the governing differential equations must be solved 
numerically. The algorithm utilized by the Patankar-Spalding computer 
code produces a solution that represents the hydrodynamic nature of 
the flow field extremely well. It has been found, however, that refine- 
ments must be made to the chemical reaction portion of the code in order 
to produce meaningful results. These refinements take the form of the 
Pratt-Wormeck chemical subroutine which has been incorporated into the 
Patankar-Spalding computer code. 

The flame shapes for pressures of 1, 5, 10, 20, and 50 atm have 
been computed and agree remarkably well with experimental data. There 
is a noticeable shape change that occurs with increasing pressure and 
believed to be a result of buoyancy effects, which become increasingly 
important at higher pressures. 

The species concentrations computed for the flow regime at hand 
seem to reflect the chemical aspect of combustion fairly well. The 
concentration profiles do not reflect much dependence on pressure, however, 
except in so far as position is concerned. Since the chemical mechanism 
incorporated into the program did not vary with pressure this is a some- 
what expected result. It would be of considerable interest to change 
the Input chemical mechanism with pressure and to note any resulting 
changes in the general shape of the concentration profiles. 
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APPENDIX 
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bEN02b86 

btN026ti8 

b£NU2690 

btN02oV2 

b-iNo2699 

btNU269fo 

bkNCXOVb 

bkN02700 

bEN027i>2 

btN0270H 

btNCe7 06 

btN0270b 

btlMG47lG 

bbN027 14 

bLNo^. / l«t 

btNG2/lt> 

btNLi:/ ifc 

bhNL474L 

btN02/44 

GEnol ? i 

beN'J4 1 e6 

bENC2726 
bcNC2730 
btNo27 32 
btN04?39 
bENC2736 
bEN0273B 
bt:NO27T0 
bENG27 94 
bENG2?99 
GE NO 27 46 
6fcN0274b 

btNG27bG 
bbNC27b2 
bEN027b4 
bENU2736 
bEN0273b 
GtN027o0 
b£NG2762 
bbNd4764 
bt NO 2766 
GEN0276O 
bfcN02ir70 
bENU27?2 
bfcN02774 
GEN04776 
btN0277t 
btN027bG 
bENG27b2 
GENt>47b4 
btN027b6 
GENU<:7bb 
btNG2/9G 
bbN04792 
ofcNL4 /Va 
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US/^bO FOKIKAN 




oiai 

0 1.03 

0 1 u3 
0 1 Go 
OiJ7 

Oloa 

UiO^r 

uiii 

0il2 

GlA*t 

Giib 

Oiia 

J X i > 


o k.- 
vX<^i. 
OXiL^ 
0*2^1 

Ol^H 
0X4:3 
012? 
0X2S 
G x3 i 
U 1 32 


1:>M 0133 
Ibrt 0x3^ 
ISrJ 0x33 
ISN 0X30 

13N 0x37 
ISN Uiia 
ISN 0X40 
XSN OIhX 
XbM Ol4x 
IbN Ul4J 


XHX-UAa:^iTK-XLul 
XtU^T A 

ih lAHi.LlxX.O) oO TU i4i 

Gutsi rxtwc 

MKX ItrlotlHGl KvAHl 
i^UKHA I ( XnUf XOX«42nnlA lUKE 

■ k*.x. aT. .-—V'Of 'k*. • a ^ k . '■r 


ItMHtKAIUKt hAlLtti lU OUNVXRGE IN 


runriA ■ I X nuf xoAfH-xinnx A I OKr icnrcKAiuKC rAxt.x.t> lu ouix 

Xi 2 tiXH lTtKA71L)N:st3Xti/Hl tnP oihFtKtNCt - * 4!:Ri)x2» 3/ | 
LUNTINoc 


XI- tLUtaUbI WKXTtlbtl42> «32iX)tK~itNtliSM,lK 

Fukmai i/3Xi2ahLum>LtTc oyraaUbriUN toT INATc/ iPiiux0.2/iPlUii0w2/ > 

XF inOubat-WaxI Ab^lON 900 ICl n£X1UK 

Xf- 1 Moot: .to >2) 300 TO Nt-XTOK 

XI- tliOut; «Gl: ».;> I Ab^lOfM 2u0 tO NrXIXJh 

AtiloiV X fk lu iVcxTNo 

o U f U ~ 0 0 




-OAKDAOi ciliMftfni (oGKUL'N ANU MCbMucI 

l^=3a^t•G 

i/LM-uoAl iNb I 
OL X7x X-XvNb 
324 X 
3M=u.l 

Xf- IMUUCatO.l) AbblON VOO TU NtXlUA 
XF XrtUuE.hU«2l AbbXUN 300 lU NfcXTOiX 
Xl^ tmjut: .Gc AbblGN 200 lU NcXiUA 
A33XGN bCG It! NbXlNG 
00 ( U :>00 


«« 

»» 4 9 9¥ 


¥♦ ¥4 ** ♦¥ JF¥ ^9 OHAPTtK 2 *9 *9 99 

¥9 ¥¥ ¥9 ¥9 ¥9 ¥9 9¥ GMARTtK 2 «9 99 ¥9 


¥¥9KXNcTXt9¥¥ 

bkOiXON XLtK NXtMfclXC solution hKOM AulAbAliL twUlLibRlUM trSllMAlbS 
ihujt: 3 AMU 4 ONLY) 

’ NfcAK-tUUlL.lbKlU«>1 SOLUTION TKlNtlXL rtlfM trtV-X*GlJ“3 aO/CU M-StLl 

^00 LcuuXL^aFALSt* 

Xx-0 

£MV = l.t>0 
XLu-cnv 

. lNLMtA:>t mNOK SPtLXLS eHOrt lOOXLXLKXOM bSiXMATcS 

OU 20 X i=i.NS 

Xh lS2tl I.LT.LMMLLI sailJ^SriALL 
20 X OONTXMUt 

assxgm 2^0 1X1 nlxtox 
assign 2X0 to NcAXNb 
bU lu :>00 

r^^lLUKc UN NtAK-cUOiL HllH trrtV«ALO* UtuKtSc CrtV 0>i AN LRUfcK 01=^ 

* — - — lout ANU AlthMPI A.bAXNt XlhKAlXNU IliX* mIaY UK TU 1 XNbS 


UAIt Ta,lVV/l0*3l«20 


bbN02796 
bbN02?Vti 
btN02aOO 
UbN02a02 
bcNu2604 
btNo2bC<> 
btNO2«0b 
bhN02ai0 
bfcN02al2 
bbN028l4 
btN02»16 
bENC2bib 
btN02fi2G 
btN02B22 
bbNG2b24 
btNt4ia2b 
UtNU2b2a 
U=N02a3U 
utNOib j2 
ucNG263*t 
b!rNi2a3t 
btNG2b3b 
bbN02b4C 
btN02b42 
bl:N02b44 
bEN02846 
bENO^ONto 
bENO2bi>0 
ULN02ad2 
btNQ2aS4 
bEN02a36 
A3ENu2b3a 
ubN02ab0 
bLN02b62 
bt:NG2Bb‘«^ 
btNOZbeb 
btN02bba 
bEN02a70 
bbN02872 
bfcNC’2b74 
btN028 76 
bENOzava 


Oi:N02BB2 
bkN02864 
bbNO^bbb 
btNU2abb 
btNOZavO 
UtNU26V2 
bLN02b94 
U t NO 2b Vb 
bbN02avb 
bt N02 VGO 
UtNU;'9tJ2 
utNX-4:sG*t 




LEtfbL 

2.2 

iStRI 

76 J 

LKcK 

ISN 

G 

c 

210 

trtV-crtV»o.i 

iSN 

Ul^> 



ALU^tMV 

ISN 

Gl^to 



1 /-IX+i 

ISn 

Qi^l 



iF liA.dw«i21 GO to 610 

ISN 

Ol^V 



TA-I SAVc 

ISN 

0130 



Oo 2il l=liNi 

ISN 

OlSl 


211 

S2i I I =SSAVt 1 1 1 

ISN 

0132 



ASSIGN 230 TO nEaTUK 

IsN 

OlSi 



ASsIgN 210 to NcXTNG 

ISN 

0l3*t 



GO lO 300 


c 




Fukiram n nXiaNuto 


C— — —rlAVc IVcAK-cWOlL bUl.(JllU!4t ^U t-lKSl iKY ulRcClLV lO 
U— Kbk*OJ.Ki.L; bULClitiU oiVtiM fcHY 

c HAV;- i>(.LAK-cUUi.l. iULUilUM, hlKliT Tk f UlKtit-lLr IlI 

L J<t;k.^lhcU 5ULoTlUt^ AI ttiV/tlM tMV 


t 


ISN 

0133 

23L. =MV = tMWSxV 





ISN 

0 130 

IF ti^OGt*rw 

»3l ASSIGN 

9u0 

iU 

NcXTcr. 

ISN 

0 i.30 

i r 1 MU U6 • i: I.J 

•SI assign 

300 

TO 

Ac A Ton 

ISN 

Jl60 

assign 230 

fU NcXlNG 




ISn 

Olol 

GG lU 300 







Ltb l ^IN 

UbT Aiiy 


iSN 

li»N 0ic>3 
iSN 0163 
1 0 1 6 o 
ISN 0l67 
ISN Oloo 
ISW 016V 
ISN 017o 

ISN on I 

ISN 0l/3 
ISN 0X73 
ISN Oi/6 


c 

C bkANCH. MaKCHING^** 

L flAVt A KXNkTiC SUUUTIUN BUI aT tHV .Li* eHVSAV. SlAKi Al 

C K.NOMN SULUTiUN ANU iNtKEASb bHV BY FaCTUR TU HUYF TDhaKOS 

i, A SUlN ThI:Kc» if SuCCbSSFULt RbFbAT UNTIL bHVSAV IS KLACnFUI 

c Nlil 6UCCLSSFUL SIAKI HALF INltRVAL SbAKLMlNG UbSCKiBtU UkLQW 

C 

2i>0 c«V=cMYAl.Lll 

it- ItNY.bl .tHySAVI brtV-FMVSAtf 
XNl=tMV 

lAsO • 

»k=ISAV£ 

UU 231 l>ifNS 

23i 1 l-SSAVtif i ) 

ASSlbN 2?0 10 NkAtOA 

It- (£HV.bt:.l:NVSAV.ANlit»KQ0E*ba.3l ASSIGN 900 TO NFXTOK 
It- (fc:f1V.Gt.UiVSAV.ANti»n(Jut:.tU»<»l ASSIGN 3G0 YU NbXlUK 
assign 2 70 I Li NLAING 
GO 10 >00 


L 

G ♦♦'►ttALh-lNl tKVAL ScaKOHINg*** 

C nAVL SULUI ILN AT AL.U OUT NOT aT AHI t H'tNLb STaKT InT^KVaL 
L StAKLniNG btf sbtTiNG bHV TO Itrfc LOGAKi TritUL AtfbKAGi: | 

L Ih XlctvAllNG NORt itiAN ILn TlnLSt URHiNAic. 

ISN ^ 

ISN 
ISN 
ISN 
ISN 
ISN 
ISN 


Oi77 

2/0 

IA=1A*1 

Ol7tt 


rX=^TSAVfe 

OlTv 


UU 2 7 1 1 =1 « NS 

0 160 

2/1 

S24 1 1 *SSAVt 1 1 1 

Gj.£il 


IF llA.Gl.iO) GU 1 

Ulbi 


c<3V-uSGK i {ALU^Xni ) 

OltiA 


AHl=tMV 


UAIfc 7ti»l9V/l0.37«20 

* ■ 

GfcN029G6 

GbN0290ti 

GtNG29T0 

GbNt,2Vi<: 

GbNOclVlA^ 

GtNuiViO 

GtrN02Vlh 

GbNc2v2U . 

GbN02922 

GtNG292A 

GtNG2'»26 

GbNi32v2b 

GtN0^9^C 

GCN02932 

GcNU2vS*« 

t»cMi'^v36 

totNOl'Y 

GBN02V‘tO 

GtN029H^<e: 

GtN029*tv 

GbNbi29A6 
GtNU^VHd 
GcNO^VSti 
GENOi’VSi 
GbN02V3t 
GbN029So 
GbN029 *>a 
1 F Gbi\l0^960 
GbN02962 
GENU2964 

GbN0296tt 
GbN0297U 
GENQ2972 . 
bbN02v7<4 
GbN029?6 
GbN0297fi 
GcN02v80 
GEN029tt2 
GEN02Vb4 
GEN029ti6 
GEN02988 
GbN02990 
GbNU2992 
GEN0299^ 

GEN02996 
GbN0299b 
GENOsOOO 
GtN03002 

GtNii^OOA 
v»FN03006 
GtNOsGOE 
GtNO^OlO 
Gbn*0sGl2 
GcNOO 1«* 
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isfcpr 

76 > OKtK OS/;160 KUkIkAN H 

OxTENUtU 



UAlt 78.199/10«37*20 

ISN 

Oias 


ASOION 2SO lU NcAim 








GbNG30i,6 

ISN 

OiSo 


ASSlON 2/0 lu NcXiNO 
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ISN 

UI87 


66 1 U Mi'O 








GtN0i020 
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bkN03022 



0*** ♦♦♦ 



OHaP I tK 

3 

♦ 9 ♦ 

*99 

Ot NO 30 29 



0 

*♦#• 



OMAK 1 OK 

3 

9VV 

9*9 

Ot;N63046 
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bEN0302a 



0 

***NCN-AutAOA 110*** 








0EN03030 



c 

StulluN ^yK N6 N*-aU1aqa 1 iC oULuT lUNS I^kQm a01a6aT1C 

cSTIMAIOS 


Ot:N03032 



0 

IMOOt 2 AND H GNLVl 








bbN03039 



0 

IKt UlKhOiOV fO LibTAlN NUN-AU lABAl 1 0 S 

UOUTiJNI IF NUT 

SOCLES SFUL 

■ t 

0tN03G3b 



0 

SlAKI mAOF'’lNrt:KVAO 6 OAL 1 NO I-kOM Tnt AOlAbATlC SGOUllON 

6V 


GEN03038 



t 

UcHIMNO A SOHLlNb KaO 1 OK FO (0 

• O'"! • 1 1 

1 U NULl IKOY 

Ih4 

NUN"AUlAbAI iC 

GtN03090 



0 

ftivK 1 m 1 IN The cNekuY cuuAlloN 

IK SPeOc* 






ObN02G92 



0 









OeN03099 

iSW 

t i da 


260 Lp,L»XAb = ,|-ALl,t:. 








6tN6309C 

I;>N 

0 AO-y 


A t.U=0 • 0 







* 

bbN030**b 

iSrt 

oi Vo 


Art j = 1.0 








6feN63CbO 

l^N 

0 1 V A 


Hx^-X.O 








OENviub^ 

ISN 

olvjj 


1 A=0 








6tNC3l; s^9 

1 sn 

ti J. VJ 


ASSloN IL NcXILa 








6CN020bfc 

ISN 

OlVn 


AX.Slui<i iSw Iti NcATNo 








OtNOiObti 

I^N 

0iV9 


m O I U Zt w 6 






‘ 


OENO ^060 



L 









0tNC30b2 

ISN 

Ga96 


320 IF IFU«cM.1«Q| OU IU VOO 








OcN03069 

ISN 

0196 


ALu-Fw 








ObN03Cob 

ISN 

0X99 


Fg^i.ti 








OEN0306b 

ISN 

0400 


AHi=l «u 








ObN030lO 

I;>N 

OcOl 


Xa=0 






- 


OEN030 74 

ISN 

O 4 Ja. 


GO IU ^10 
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0EN03G 19 



0 









b£NG307b 

ISN 

040:1 


330 IX=1X'*1 






4 


bEN03b7t 

ISN 

020^ 


Ih tiX.OlaXOj OU To 2AO 








btN030bO 

ISn 

0400 


TKs=T:>av6 








O|;N030tt2 

ISN 

040 7 


UU i'liNS 








OtN63Ub9 

ISN 

O^Oo 


6i3x S 4 U )=SSA9t 1 1 * 








bEN030bG 

ISN 

0409 


rv=0 *3* i XOU^Xni I 








bbN030bb 

ISN 

04io 


Xril^'Fw 








btN03GSQ 

ISN 

0211 


GO IU :>XC 








btNO3092 
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ObN030V9 

ISN 

04i4 


296 OUNl iNUt 








bENGS096 

ISn 

021^ 


FO=i* 








bbN03GVb 

ISN 

02l*t 


6U TU c30 








bo NO 00 



c 









bON03l02 



O'#-*-*-*- ♦■♦♦V ♦♦♦♦ 


CHAP 1 tk 

9 

»«♦« 

***OcN031C9 



0***# ««*«■ *«*¥ ««** 


CHAPTeK 

9 

9999 

9«9btN03106 



0 









bEN03lLb 



0 

***FAi6UKt tAliS*** 








bfcNG3ilO 



U 

rAiOcO tUulL UK KlNtllC 6U0N UK 

tout V 

KA 1 J U 

UUlSIUc 

1 0. 

itlOl 


bEN02ll2 



C 

KtlUKN AulAbAIXOt NON-KcAOTtO NlXlUKc 

KKuPcHl 1 cS 





btNu3il9 



c 









bfcN03ll6 

ISN 

021r^ 


960 Srt-U.O 








btN031 ib 

ISN 

0216 


UU 90 i 1-1* NS 








btN0312G 

I:>N 

02x1 


S2ii I=6i (i) 








btNOil24 

ISN 

0218 


SM = SN+62 i I 1 








bcN0jil«.H 


US/:»6D hUKTKAN t% tXlkNUtrU 


llATt TB*iV9/iO»3/*20 


Q^ilV 

0 «: 4 i 

0222 

u 2 iJ 

02 *;^ 

Ol£a 

0 ^ 2 e» 

022 ? 

0220 

022 v 

023 o 

U 22 i 

U 4 .:>^ 

G 2 ;J^ 

0 <i 3 <^ 

0 *l 37 

02 


ISN 02«»x 
ISn 

02^4 
l:»N 0^40 
Ion 0247 
ISN 024b 
ISN 02 SO 
ISN 02SX 


ISN 0252 
iSN 02S4 
ISN 02o:> 


CUNi XNUi: 

rN=-iooo«u 

A Lit- IK 
InL KO*i 
7 k 1 NV=x .OU-j 
UL 400 l=l»At. 

Call mlhs 
n bUH— u • u 
OU 40*^ K=-A»No 
hSOW-nsOrt+rtulK >»S2iK* 

CliN I INlc 

T K=: I K.^1 i .C^IhSOao»KGASlN^l NlNV-HSOrt l/CHSUrtI 
Tr^lNV-i * 0 /lA 
Xni^UALOl I^~ aLU» 

X lU— I X 

ir iXlit.Ll.X.Cl GO iLI tO-* 
uLNliJVoc 

WKilblLfXnCl ifAHI 

Clis ix nlc. 

isr:Lie=r'#^/ IKOAO* iK'f'SMI 
OC (u VOL 


o 

c 










CfiAI'T ch 
OrtAHTfeR 




GtN03l26 
beN 03 i 2 U 
otNOSiii) 
GENUOI32 
OtN 03 l 3 H 
OtrNOoliO 
GtNOOlSb 
Ut;N03i40 
btN03l42 
UEN03144 
btN03l46 
btNOS I 4b 
b 1 i:N 03 l 5 t 
bbN 0 O 152 
Gl:NC3Xd4 
OtNOoiSo 
OcNC315i> 
o^:Nuol^0 
ubNOOl 
OtNOolbH 
bbNOolOb 
bbNOoxot 
^0fcN03l ?C 
«bbN031Y2 
bEN03174 


»*.*r>KtibLbrt UbLL CALL ¥0 sPbLt**^ bbNOBlTb 
lAKb tnc bSIlNAtES bbNERAlbU IN ChAVlEftS lf 2 ANO AIIENRI A SOLUllON GbN 03 I 7 b 
vUTn t-OLL bOUATiUNS. II- SuCCbSSEOL* uKOATt TMk SAVb ANSNbHbS Wlln InfcGbNO^jlbG 


OUlUIION mNu KcIUKN 1u blATfcKcNI NOHltbK NbXTUlL* IF NUl t thb aCI lUN GbN03Ib2 
UbRclMub UN MHtrlhtK AN bUOlLbH UK KiNbTiL SULN IS bOUbnT. rMibbU bbN0OIb4 

bOOlL bULNi KblUKN ¥U SlAlbHbNT NUliBEK NEXlNb* MhlLb FAILURE IN A GbNOBIbb 
MUcIlL oOLN HILL Be HjJLL0«4bU BY AN Al'TbMPf WITH LNRb-|- — - — » kT*^0.0 bbN03l8B 
ANu bAne LBllMAIeS. bblllNb K1=U»W IriPLitb IriAl a UHANGE IN lEiiP GEN031V0 

FibLU NAS NO ebi-eCT UN SPtrClbS lilSTKItSUT lUN hUK ThaT PAKTICULAK GfcN03IV2 

IIcRaIIUNi but uObS ALLOW Tub SPcLibS CnANbbS fU iNpLubNCb iNb fbNP bbN03iS4 

UhANbb . HaKIIAl UbLUOPLlNb OF Tht bNtKbT EWUAUUN* IF SUCCfeSSFULGEN 03 i 9 o 

riiln LNKb=r t kcPlaI WITn lNku^^I IFULL bOUAnONSH IF bTIlL NU GUUU* GbNOBlVb 
K::IUKN tU blAIbNbNI nUMBcK NEaING* bbN03200 

GCN03202 
GbN0320t 
bbN032C;6 

•U bbN0320b 

GFN03210 
bbN03212 

70t2I0t2B0«27Ct330«60O| GENGb2i4 

bcNU32|ib 
GbN032ld 
bENG322C 

IbKNbt StNU TnIS bbN03222 

I LALCOLaIIUn GbNUB224 

GbNG32?o 
btN0322b 
btN03230 
bEN032 32 

UPl Art OAVl ANbWLK::^ AfML* LUi'<T INUL* bbNoji23<*^ 


Call SPcul 

IF iLLUNVbl bU TU ::40 
SULUIiUN FAIlEU try Kl-O.U 
IF (LNKo) bO TU :>^0 
LNKb- . I KUb • 

GU 1U NcXrNbt liO0tI30*i70t2IGf2B0«27Ct330«60O| 

If tUbuUlLI bU lU MU 

LNKb-»FALSt* 

Gu iU :>0'0 

MAVb ioLN uUl WIlNUUr kf ItHMbt StNU TnIS 
>UlN AO eGllNAleS wllh Kl LALCOLaIIUn 
If iLNKbl GU To 550 

LNRb^a IKOc* 
bU tu :>0G 


0 — :>ULur ll-N 


:>ULLlo 5 F JL t 


■oQ 


09 

Ir 

3s 




LEVfcL 


IStRT /o) CKtK US/3&0 EUkTRAN H tiXIcNOtU DAfb 76 .199/10 




C 

btNG323e> 

ii>N 

O^DO 


3SC TSAVE-XR 

bcN0323b 

XSN 



Uti 3BO 

bfcNC32At 

ISn 



3oC c>^AVt 1 i J i 1 

btN032-t2 




C 

btNbJ 2 **<* 

i:>N 

SV 


bU ro ucAlURf {200f2iOf23U*3CC »32G*VUG| 

btNCj 2 AC 




c 

bbN032Ab 




^t***** CHAHI’EK 0 

i>bN032d0 




c-***^^*^ 4^****^ CHARTER 6 

bENC32S^ 




t 

bbN032iv 




C v»2ckKuK rtcibAGtS*»» 

btNCb^ &<> 




0 

bbN032bt 

1 S(M 

0 .>o 0 


Bltb WRX t c 1 ti t i 1 

btN032c0 

lin'i 



oCj. *3uH hAlLUKt: Tu biNu trOUXt S0uN..«i 

bcNu32t>2 




A2tnAtfl> xNt_Li RKURS KtiLiKNtu/l 

obN032 tjH 

i SN 

0^ Di:. 


bU ID ts3o 

bEN032oB 

I^N 

0 ]clC) ^ 


Ox>^ R.KiltlCfCXXl 

bbNC-32C)b 

ISN 

U £.0 <v 


Oil rCKMAfl lhu»lu>.»3l*fH^»»^l * 33 H EA1LU»-i- TU EJRu NtAK-cOUlL bOLN.*.* 

bcND327L 




A2*»tiAWt> liVLcl pRLiRS Kc 1 Ot<iMC 4 ,>/ 1 

OCN03272 

i:>N 

Gko3 


0 1 . 1 ij 0 3 C 

btNu32 Im 

l5iM 

Oitic* 


02 u wr 1 

bbNl*S27c> 

ISn 

t 


o*:l pUcwMAI iinv#i»jAtJlHn»>fr^*l»bih FAlLuKt lU UtifAlN RlNcTlb bULN ARTcK 

fbbN0327b 




A-rcn lc:N IN 1 ct\ VML rAL ViNis. • .A Vb XNLrl RMjRb KEIURNcU/I 

bbN032b0 

iSN 



bU to B3G 

bbNC32b2 

ItiN 

O^toV 


bbC WKllElBtbbll 
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